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STRESSES DUE TO THE PRESSURE OF ONE 
ELASTIC SOLID UPON ANOTHER 


I. InrTRopucTION 


1. Practical and Historical Aspects.—The question as to what 
stresses are produced within the area of contact when one elastic solid 
is pressed against another has for a long time been recognized as one 
of great importance. That a solution of the question has a definite 
practical application may be realized when it is considered that the 
calculation of the stresses encountered in ball and roller bearings, 
bridge rollers, rolling lift bridges and wheels on rails is impossible by 
the common formulas of mechanics. The design of such structural 
parts is based largely on results of tests, or trial and error methods. 

This problem was attacked mathematically by H. Hertz who in 
1881 published a paper* “On the Contact of Elastic Solids,” and in 
the following year an extension of the same analysis under the title 
“On the Contact of Rigid Elastic Solids and on Hardness.’’ Previous 
to Hertz, Winkler and Grashof had both attempted mathematical 
solutions, but their results were either approximate or involved some 
unknown empirical factor. Hertz’s analysis has long been considered 
the basic method of solution and probably will give accurate results, 
provided the radius of curvature is large in comparison with the size 
of the area of contact. However, Hertz’s computations deal defi- 
nitely only with surface stresses, and he fails to recognize the great 
importance of components of stress at right angles to the direction 
of pressure. 

Even though by the use of a mathematical method of solution the 
various stress components (principal stresses) at points within the 
area of contact are calculated, still the interpretation of the signifi- 
cance of the magnitude of the stress components requires further 
explanation. 


2. Theories of Failure.—When at a point in a body there exist 
stress components in more than one direction, the conditions under 
which failure will occur have been the subject of numerous tests and 
speculations, and as a consequence there have been proposed four 
principal theories of failure.t By failure, in this connection, is meant 


*H. Hertz, Gesammelte Werke, Vol. 1 (Leipzig) 1895. English Translation in ‘‘Miscellaneous 


Papers,’ H. Hertz, 1896. ; 
y +For a more detailed, yet brief, discussion of the theories of strength see Timoshenko and Lessells, 


“Applied Elasticity,’ first edition, 1925, Chapter XVII. 
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the occurrence of a plastic yielding or flow of the material in the case 
of ductile materials or rupture in the case of brittle materials. The 
four theories mentioned are: 


Maximum Stress Theory 


Under the maximum stress (Rankine) theory it is assumed that 
failure will occur whenever the maximum stress component equals or 
exceeds the yield point of the material for that kind of stress, without 
regard to the effect of other stress components that may be present. 
This is the common theory and is sufficiently accurate if the stress is 
almost entirely in one direction. 


The Maximum Strain Theory 


By the maximum strain (Saint Venant) theory it is assumed that 
the maximum strain (deformation) due to the stresses acting at a 
point is the criterion of failure. This theory also has been found to 
be limited in correctness. 


The Maximum Shear Theory 


By the maximum shear (Coulomb) theory it is assumed that 
failure will occur when the shearing stress on some plane (usually 
oblique) due to the stress components is equal to the shearing yield 
point of the material. It is generally agreed that this theory is satis- 
factory for ductile material having approximately equal yield points 
in tension and compression. 


Mohr Theory 


The Mohr theory is an extension of the maximum shear theory 


which makes it possible to deal with such materials as cast iron where _ 


the limits in compression and tension are widely different. It be- 


comes the equivalent of the maximum shear theory when the tensile 


and compressive limits are equal. 


It is not considered necessary to discuss Haigh’s Strain- Energy — | 


theory, which is more recent than those mentioned. In all of the 
calculations and discussion which are to follow the maximum shear 
theory will be taken as the criterion of strength. 

In text books on mechanics of materials* will be found the proofll 


that at any point in a body subjected to stress components acting at 


right angles to each other the maximum shearing stress is numerically 


equal to one-half of the algebraic difference between the maximum — 


and minimum stress components, and that the plane of the maximum 
shearing stress bisects the angle between the stress components. 


*For example: Fred B. Seely, ‘‘Resistance of Materials,’ Wiley & Sons, 1925, page 53. 
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3. Division of Work.—In the preparation of the material for this 
bulletin, the work of the two authors has been in large measure dis- 
tinctly independent, although each has critically read the other’s 
manuscript. To V. A. Hoersch belongs the credit for the whole of 
Chapter II, “Mathematical Derivation,” which is the mathematical 
work upon which the bulletin is based. 

The laboratory work involved in the verification of the mathemat- 
ical analysis was under the direction of H. R. Thomas, who also has 
written all except Chapter II. The stress-component factors for 
Tables 1 and 2 were calculated by both authors. 


4, Acknowledgment.—This study was made possible by funds pro- 
vided by the Urinirizrs ResearcH Commission, Wm. L. Aszort, 
Chairman, for an investigation of transverse fissures in rails now 
under way in the Fatigue of Metals Laboratory at the University of 
Illinois. The following members of the Commission were appointed 
as an advisory committee for this investigation: 

G. F. Tesperts, Engineer of Structures, Chicago North Shore 
and Milwaukee R. R. Co., Chairman. 

A. Herz, Engineer of Tests, Public Service Company. 

J.S. Hyatt, Engineering Assistant to Vice President, Chicago 
Rapid Transit Company. 

A number of meetings have been held with the advisory com- 
mittee during the course of this work and they have given much 
helpful advice. 

The investigation has been carried on as a part of the work of the 
Engineering Experiment Station at the University of Illinois and has 
been under the general administrative direction of Dran M. §S. 
Ketcuvum, director of the Engineering Experiment Station, and of 
Pror. M. L. Encer, head of the Department of Theoretical and 
Applied Mechanics. 

All work in the Fatigue of Metals Laboratory is under the direc- 
tion and supervision of Prof. H. F. Moors, and to him the authors 
wish to acknowledge their indebtedness for many helpful suggestions 
and for sympathetic encouragement during all stages of this study. 
The careful and painstaking work of Mr. N. H. Roy, Research 
Graduate Assistant, was of great value in the laboratory work in- 
volved in this problem. 

The idea that the shearing stress in the interior of bodies in con- 
tact might be a critical factor was suggested some years ago by Mr. 
T. M. Jasper, who was at that time Engineer of Tests in the Fatigue 
of Metals Laboratory, University of Illinois. 
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Il. MatruEeMaticaL DERIVATION 


5. Introduction to Theory.—Two bodies are considered to be orig- 
inally in geometrical contact at a single point and then pressed to- 
gether so that there is contact over a small area called the compressed 
area, bounded by a curve called the curve of compression. The theory 
of the stress set up in either body was first developed by H. Hertz.* 
Hertz completely solved the problem but left his solution in terms of a 
certain Newtonian potential function. The integral involved in this 
function was not expressed in terms of standard elliptic integrals. 
The innovation in the present work{ consists in doing this for points 
directly under the center of the compressed area. A rather condensed 
account of Hertz’s theory appears in Love, ‘“‘Mathematical Theory of 
Elasticity,’ 3rd edition, 1920, pp. 190-196. All references to Love 
are to this book. The notation of Love will here be followed. The 
author will present what he hopes is a simplified account of the Hertz 
theory, followed by his own elaboration of it. 

Consider the dimensions of the compressed area as small in com-' 
parison with the radii of curvature of bodies 1 and 2. Body 1 may 
then be considered in the unstressed state as an infinite solid 
occupying all of space on the positive side of a plane z; = 0. When 
body 2 is pressed against body 1, the stress in body 1 may be regarded 
as due to a pressure P’ distributed over the compressed area regarded 
as an area in the plane z; = 0. One system of coérdinates 21y121 will 
be used for body 1 and another system wey2z2 for body 2. Before the 
bodies are pressed together, the x, and 2, axes coincide, the y; and y2 
axes coincide, and the z; and z2 axes are directed toward the interiors 
of bodies 1 and 2, respectively. When the bodies are pressed together, 
the surface, z: = 0, is deformed into a new surface, 21 = wi, where w; 
is a function of xy:. In general, a point x1y121 is displaced to a1 + wi, 
Yi + 01, 21 + wi, where ww, are zero when the point 21y12: is at an 
infinite distance from the compressed area. For body 2, a point xoyoz_ 
is displaced to ao + ws, yo + v2, 22 + we in a similar fashion. The 
original point of contact of bodies 1 and 2 is the origin (0, 0, 0) in both 
systems of codrdinates. When the bodies are pressed together, this 
point is displaced to 0, 0, w: in the one system and 0, 0, we in the 


*An English translation of his theory may be found in “Miscellaneous P 1 
(Jones & Schott, London, 1896, p. 146). ee aires tak hg 

{Since this was written there have come to the author’s attention two bulletins written by Prof. 
N. M. Belajef, Professor, Mechanical Laboratory, Institute of Ways of Communications, Leningrad, 
Russia. In these bulletins Prof. Belajef solves this same problem, also making use of elliptic integrals 
for this purpose. Although expressed in different form, his equations for X. x Yy,and Z, agree exactly 
with equations (59). However, Prof. Belajef apparently considers that a tensile stress within the 
body is numerically equal to the difference in stress components; in this bulletin the shearing stress, 
equal to one half of the stress difference, 1s considered to be the critical stress. It is interesting to 
sae shat 7 makes use of the mathematical analysis for computing the stresses in a railroad rail under 
a wheel load. 
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other system and the displaced point is the center of the compressed 
area. Then a = w; + w» is the distance by which the two bodies 
approach each other when they are pressed together. 


6. Definitions.—The following notation will be used: 
Ai, #1; Av, we Lamé elastic constants for bodies 1 and 2; yw, we are the 
rigidities (moduli of elasticity in shear) and \, + 24u1, 
he + 282 the moduli of compression. 


E,, FE, Young’s moduli (moduli of elasticity) for bodies 1 and 2. 

O1, 02 Poisson’s ratios for bodies 1 and 2. 

p1, p's Reciprocals of the principal radii of curvature of the 
surface of body 1 at the origin. 

po, p's Ditto for body 2. 

&1, m Distances from a point x:y121 to the planes of curvatures 
p’, and p; respectively. 

£, Ne Distances from a point %¥222 to the planes of the 
curvatures p’» and pz respectively. 

€ Angle which the planes of the curvatures p; and p2 
make with each other. 

LENG Normal components of the stress across planes perpen- 


dicular to the z, y, and z directions through the point. 
Tension is positive and compression negative. 


ee Z, for points in the compressed area. 
iP Total force with which one body presses on the other. 
K Ou; Ov; OW, ‘A U2 Ove OW2 
Pe a aL ae ae ay AG Se os 
0” 0? 0? 
2— 
M Oe ie oy? ss 02? 


7. Action of a Force Concentrated at a Point.—Consider first the 
simplified but ideal case of a force P concentrated at the origin instead 
of being distributed over a compressed area. The strain at a point 
x, y, z in the solid is specified by the following equations written in 
different form but equivalent to those of Love.* 


Pz 2 (2) P i. 
eg dan dg r) 4n(A + p) Ox SSG 
PZ wae Kl: Ie ) 
ee ry. eho ) 
: Amp oy (*) 4An(X + pw) dy loge) 


w= log (2 + r) 


EZ no (*) P(X + 2u) @ 
Amp dz Amu(X + p) 02 


*Love, p. 189, Eq. 35. 


If 
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where raVoety+e2. 
Then 
Pz 1 i Pad () 
Jo ba ol ee Se eee See ee eee} | 
4a ye (*) 4m(X + p) pay a 4nrpdz\r te 
P(A + 3n) 0? 
Amu(d + py) 02? Be 2) 
Since 
ie sipcenteattis i 
V 7 = 0, Vlog (¢ + 7r) = 0, an oe og (z+ 1r) = Es 
this becomes 
183 @) ff il 
ee eae td (2) 


The stress may be determined from the strain.* 


x eae vA Sie 7 NA jae 3) 
i SONA Bay? ae sitar Tir See + 2u ( | 


| 
No Etc al es | 
Using (3) and (4) we have: | 


r 
= 20; s 
x oP a @ REO" (*) (1 — 20)P 0? , 
= mg O2\r) Qe aa? \r] Qr Ox? SVN) 
y 22 ()_ Pa) _Gaa0r 2 6) 
ve de\r) — Qe dy? \r] Qa oy? log ot) | 


pan dah Oa or 8 
"Qn d2\r) — Qa 022 \r 
8. Solution of General Problem in Terms of P.—Let us return to — 
the actual case of a compressed area. The effect of P’ distributed — 


| 
over this area is the integrated effect of forces P’dx’dy’ over the com- — 
pressed area. The components of displacement are then obtained | 


from (1) by letting r equal V(x — a’)? + (y—y’')? + 22 and taking 


*Love, p. 124, Eq. 8; p. 101, Eq. 21. 
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ig 


Yolies Of Stress CHF 


A0 OF OF GE QB V2) WE 14 16 LO EO 
Yatwes of 


Fig. 1. VARIATION IN STRESS CoMPONENTS IN DISTANCE FROM CONTACT SURFACE 


the integral dz’dy’ of each term over the compressed area. Introduc- 
ing the symbols 


P’ 
o= [fF dz'dy’ and x = ff log(z + r)dx'dy’ (6) 


and interchanging the order of differentiation and integration 


z@ 0d 1 Ox 
Oe Anu 0x An(X\ + pw) Ox 
pane z Od a 1 Ox (7) 
Any dy = 4a(X+ #) OY 
z dd A+ 2u Ox 
et 4mru dz Anu(A + p) 02 
By similar process from (2) and (5) 
1 Ke) 
A (8) 


~ Qe(X + pw) d2 
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e « Oo z 0d Teo 0? x 
ee nw 02 20 0x2? Br Ox? 
0 20 1 — 20 0? 
ype fue one * (9) 
w 02 27 Oy? 2x oy 
1 0¢ z 0h 
Lo 


Qn dz Qa Oz? 


d 
Since es log (e + r) = (1/r) and since X, and Y, must evidently 


vanish at infinity, we may write: 


02x je 0? x {2 
Gk sess ee 10 
Ox? Ay rice o* Oy? OU ae ey 


It is clear that Equations (9) contain the solution of the problem 
when the eurve of compression and the pressure P’ within it are 
known. The bodies have heretofore been regarded as having plane 
surfaces. To determine the curve of compression and the value of P’ 
as a function of «’, y’ the surfaces must be regarded as curved, and 
afterwards they will again be regarded as plane. In other words, the 
curvature of the bodies is regarded as having no other effect on the 
problem than the determination of the curve of compression and the 
value of P’ as a function 2’, y’ within this curve. 


9. Determination of P’ and ¢.—The equation of the surface of 
body 1 in the neighborhood of the point of contact is 2: = 1% pié?2 + 
VY p'1n:” and similarly for body 2, 22 = % pof.? + % pons. Then by 
the formula for rotation of axes through an angle e: 


& = £ cose + m1 Sine, no = —£, Sine + 71 COSE 
Hence 2=2,+ 2 = A’é? + 2A’bin1 + Bn? (11) 
where 2A’ = pi + po coste + p’o sin’ 
2H’ = (po — p’2) sine cose (12) 


2B’ = p's + po sin2e + p's cose 


By a rotation of the proper angle from the £1, m1 axes to the x, y axes 
we get (Appendix A) 


2 = 2, + 2, = Ax? + By? (B 2 A) (13) 
where A, B are the roots of the equation in x 


(2% — pi — p2)(2% — p's — p’2) + (p1 — p'1)(p2 — p's) sin?te = 0 (14) 
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If e = 0, then A = 14(p1 + 2), B = 14(p's + p's) (15) 


When the bodies are pressed together, let the points (x, y, 21) of 
body 1 and (a, y, 22) of body 2 come into contact within the com- 


pressed area. Then 2; + wi = —(z. + w.) + a. The minus sign is 
due to the opposite directions of the z; and z: axes. Hence by (13) 
Wi + We = a — (24; + 2) = a — Az? — By? (16) 

For points outside the compressed area 
Wi + wW2>a— Ax? — By? (17) 


For a point in the plane z = 0, we have from the third equation of (7) 


At + 2u1 € | 
= 18 
st Amui(rr + w1) \02/, =0 i 


i 7 ie 
2 es oe : Tear , Dah? 1 
Be s/f Ae log (2 + r)da’dy P'(1/r)dz'dy (19) 


It is known that* 


(3r1 a5 2441) M1 
Io eee SORA cae 20 
ee Nee i er ee SOE AL) au) 
and therefore 
a 2 
] O1 ae Ay + Qu (21) 
Ey 4ui(Ar + 1) 


Using (19) and (21) in (18) and letting ¢ = ¢. when z = 0, we have: 


1 —) o1- 
Wi = 0i¢0 where 3; = aig 
(22) 
il an To" 
Similarly We = Jed. where 32 = mii, 
Substituting w; and we from (22) into (16) and solving for ¢. 
1 
= — Ax? — By? 23 
ange OE) (23) 


*Love, p. 124, Eq. 8; p. 101, Eq. 21. 
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Judging from (13) it seems probable that the curve of compression 
is an ellipse. Now, as defined by (6), ¢ is a Newtonian potential 
function due to matter distributed over the compressed area with 
surface density P’. Hertz recalled that the potential at a point 
z, y, 2 in the interior of a homogeneous ellipsoid is of the form* 
a — Ax? — By? — Cz? where the 2, y, and z codrdinate planes are the 
planes of symmetry of the ellipsoid. For the plane z = 0 this yields 
the form required by (23). Since ¢, defined by (6), is the potential 
due to the attraction of a surface layer of mass, the potential due to 
an infinitely flattened homogeneous ellipsoid must be taken for ¢; 
that is, the potential due to a mass density 


BP \ x"? fe 
BOP Pe age (24) 
gl? "2 
distributed over an ellipse = + 7) 


/ 


= 1 in the plane z = 0. 


From the known potential of a homogeneous ellipsoid} there is readily 
found, for the elliptical layer just described, 


Pale ‘be v7 y? 4 dy 
“4d, Vo ety Pty vi@+oer oy © 


where y is the positive root of the equation 


ye ie 22 
~a@ty By 7 2 8) 


f(z, Y; z, 7) = 1 


Since y = 2 = 0 within the compressed area 


d a os be a y? y 
eek ety B+tv/l@+Hnetrpye 2 


Since the values of ¢o in (27) and (23) are equal, the coefficients of 1 


x*, and y* must be equal and we have: 


7 GON ears Du oo 


We” p. 194; Peirce, Newtonian Potential Function, 3rd ed., 1902, pp. 117-126. 
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Yee, = 0 st 0) f- (a? Je z a aw 
(29) 
Be as Oo O she he 
= 4 ( iets 2) é (a? + y)” (b? + y)*y” 


Equations (29) may be solved for a and b. Hertz solves (29) for a in 
terms of one transcendental function of B/A and for b in terms of 
another such function. Tables of values are computed for these 
functions and elliptic integrals are given a mere mention by Hertz as 
the basis for computation. The problem of the solution of Equation 
(29) for a and 6 will be taken up later. 


10. Stress on Z Axis in Terms of Elliptic Integrals.—The value of ¢ 
from (25) can now be substituted in (9) and (10). Equations (7) and 
(8) will only be applied for points (2 = y = 0) on the z axis. In this 
case we have from (26) y = 2%. To find derivatives of ¢ the rules 
must be used 
Co) Oy 6 ) ) Oy 6 6 to) Oy 6 


je > A eR es ae 


(30) 


where 6 indicates differentiation when ¥ is regarded as an independent 
variable; that is, when no account is taken of (26). We thus find for 
x=y=0 
a6 Pe dy 
RI SB) CR) 


(31) 


Changing the variable of integration to w where y = a?w* and writing 


ke=-b/a. and <= 2/4 
d 
ag ” ea w (32) 


dz a’ w?(1 + w?)” (k? + w?)? 
Similarly 
02h 3P 2 dw 
ant Oa J, A Fo) (Ke? + wd Ge 
a’ BP 
ay? = we Ge ae rey oe) 


0° ef dw 
d22 a8 Je w®(1 + wo?) (k? + 2)” 
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T @rva ae (k? + £) ee 

Differentiating (33) and (34) with respect to z = af 
a = ; aan ak age 4 ¢28 (36) 
BO Leen? ou on 


dzdy2 a Otay? i at(1 + £2)” (k? + 2) 
From (10) and integration by parts using (36) and he 


0*x ” 0° 0° ee 
Pomme ei ee 


Ox? eta eS ee 


3P¢ dw ni cdg 
—~e@ J), Gta toy® Tel, C+ O84 oF 


aie 


3Pt me A rire 
eh, Cte @+o7*® — 2a —-wVI+e 
3 
+ od — 


02x ~O2q a2 nae dh ne 
Boni hae 21 a” Gray! Ge este 


taf aes dg 
t dzd0y? 


RES af 3P Hd edt 
r (1+ a is + a7) T a +) (+E 


arf ate + x2 
aot ve ae Ca ne i ac —wD Vet e 


ia =F 


Substituting derivatives of ¢ and x obtained in this paragraph 
into (9) the normal components of stress X,, Y,, Z, are obtained. 
By consideration of symmetry for points on the z axis, the correspond- 


(38) 


(39) 
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ing tangential stress components must vanish; so that X,, Y, and Z, 
are the principal components of stress. The resulting formulas for 
the principal components of stress are not written out here since the 
integrals involved will be expressed in terms of elliptic integrals and 
the equations for the stress will then reduce to Equations (59). 

The author’s innovation consists in expressing the integrals ap- 
pearing in the last paragraph in terms of F'(¢, k’) and E(g, k’), the 
first and second elliptic integrals, for which tables have been com- 
puted. These integrals are defined by 


° dé ag ee 
/ =, SSS —=B Vs eT ey nai 
F(@, k’) =f Ni awiant | BCe,*) =f 1 —k’?sin?6 do (40) 


and the corresponding definite integrals are obtained: 


eee f Lee 
Ke) = FG.) = > Ni—wtsinta ? 


(41) 
E(k’) = E (Zn) = [? \1 — k’? sin 20 do 
0 
Changing the variable of integration to w = cot 6 

° (kh? + w?)”* 
Fg, k’) =f Giz “sey wo) ;H(¢,k’) = k a4 nnd (42) 

*(k? = w?)” 
K(k’) -{ (1 at “ar os w?)” ; E(k’) = ca F (1 BBs 2) 8 dw (43) 
where k? = 1 — k’? and ¢ = cot ¢ (44) 


For the integrals appearing in the preceding paragraph the following 
results were obtained: 


sli dw ie E(¢, k’) (45) 
: w?(1 ae w?)” (k? +. w?)” A ke k? 


dw 
(1 + 0)” (k® + 0%)" 


shee a+ ae (e+ ant = = (a) ae X) | 


T Mh 
~ peng * G7 


= va Ee k’) — E(@, X) | (46) 


(47) 
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72 2 
where T= y' ane | (48) 


The direct derivation of these results is rather tedious but the inte- 
gration may be verified by differentiation, since it is known from 
(42) that 

dF (9, k’) 1 dE(¢,k’) (k2 + 52) 


de Git EP 


(49) 


The problem of solving Equations (29) for a and b may now be 
resumed. Changing the variable of integration in (29) to w where 
y = a’w’, the integrals become 2/a* times the integrals of (46) and 
(47), respectively, when the lower limit of these integrals is ¢ = 0. 


Hence 
Cees 
= ——~oaapn K(k’) — E(k’) (50) 
3P(8, + do) F 1 
B= oF le E(k’) — Ke) | (51) 
Dividing (51) by (50) 
B (1/k*) E(k’) — K(k’) 
A~ K(k) — E®) 2) 


From (52) a graph may be plotted connecting k and B/A, points on 
this graph being found by assigning values to k and computing cor- 
responding values of B/A. (See Fig. 2.) Since by (14) B and A are 
functions of p1, p’1, px, p’2 and e, this graph determines k in terms of 
the geometrical configuration of the two bodies. Adding (50) and (51) 


Bye (dy —- 0) E(k’) 
aie 27a? ke (53) 


27a* 7 
3P and using (22) we have: 


Solving for 


27a a AE(k’) AP 2rk? 
Spe ke? or a3 = 3E(k’) (54) 


where 


hee 


(3; + 32) ss i 1 — oi? 1- o2” 
A+R 2> Cee eee (55) 
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alles 
as 


Ss 
8 


S 
& 
loties OF Fe ord 


i BZ SH Ff coos a Ye 8 
VOWVUIES OF a 


Fic. 2. VARIATION IN Certarn Factors with Cuances 1n B/A 


By the theory of equations A + B must equal —14 times the coeffi- 
cient of x in (14) when this equation is expanded, i.e., 


A+ B= (01 + p's + po + p’2) (56) 
Substituting this in (55) 


1 — o;? a ; , , 
Kh = @ E, -+- E, + (p1 + p’1 + po + p's) (57) 
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When k has been found as just described, the substitution of this 
value of k in (54) yields a; and thus the semi-axes, a and b ( = ka) of 
the ellipse of pressure are determined. 

With the value of A from (54) we have: 


NOE 2k? xa? AY DA me" AZ, | 2k? a? 
@ BUR) 3P eae Ee a ee re 


Z. (58) 


Let X,, Y,, Z. as given by (9) be expressed in terms of elliptic inte- 
grals by means of the derivatives of ¢ and x in the first paragraph of 
this section and the integral evaluations (45), (46), and (47). Upon 
substitution of the values of X,, Y,, and Z, thus found in (58) 


Went 
= M(Q, a a ») 
a 
AVS 
~ = M(Q, + 0f,) (59) 
ei Goad) 
Gay OE oe 
2k? 
where M= KE (h) (60) 
1S Fe 
and OQ, = — 5 + SIF e ks EC e, kD] 
= — +14 MO/e Ee, HY) — Fle, bY] 
(61) 
1 il T : ; ; 
% = 97+ 9 7 I t SIC/h) Ee, k’) — FC¢,k’)] 


MY, = —-1+74+ S[Fl¢, k’) — E(¢,k’)] 


By these formulas, the principal components of stress can be found at 
any depth z = a¢ below the surface, since tables for F(¢, k’) and 
E(¢, k’) are available.* When tables are computed for the stress 
components at various depths (see Tables 1 and 2) it is found that for 
a wide range in values of k, the maximum difference of stress com- 
ponents, Z, — Y,, occurs at a point below the surface of contact and 


*Hancock, Elliptic Integrals; Jahnke & Emde, Funktionentafeln. 
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not at the surface. This is discussed in more detail elsewhere. It 
only remains to examine certain special cases. 


11. Special Case, k = 1.—For k = 1, b = a and the compressed 
area is a circle. The integral in (45) becomes, when ¢ = cot 4g, 


a dw ¢ 
F RegOvES ap hed Ms mad cae tele (62) 


The integrals in (46) and (47) become identical and equal to 
fs J sint ode = -¥6 si +34 ¢ (63) 
© esas), sin? ypdg = 6 sin 9 cos 49 
Letting Y = a?w?, Equations (29) become 


3P ae dw "i 30P Bi 
Ava p= 2a3 (1 + de) F (1 ab w?)? ae 8a3 (Oy + 32) ( ) 


: /_>- tes 1 choy 
since by (63) eco = [—Msin pcos p+ Yo], = | 


Tv 


Since by (41) H(0) = fp? dé = 5 , we have on substitution of this 
0 


value in (58) 


INO. 4a? AY 4a? AZ, 4a? 
a see oes 5 = 3p 2 (65) 


By the same process described for getting (59) except that the integral 
evaluations no longer involve elliptic functions, we get 

NX AY 2 : 

= = —* = = {(1 + o)(-1 + eeoty) + sin? o} (66) 


AZ, 
a 


2 
= — —sin’¢ (67) 
us 


M. T. Huber* has found the stress in this case (&k = 1) for a point 
z, y, 2 in a general position. His solution when specialized for the 
z axis agrees with the above. 


*Ann, d. Phys. 1904, Bd. 14, p. 153. 
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12. Special Case, k = 0.—Next consider the case of k = 0. 
Equation (54) shows that a = © since A and 


E(1) = f Ni sint odo =1 
0 


are both finite. A simple illustration of this case is that of a cylinder 
z, = By? resting on a plane 2, = 0. Before the bodies are pressed 
together, there is contact along the x axis. When pressed together, 
the two bodies are in contact over the strip bounded by the lines 
y =band y = —b. Let P be the integrated pressure across the 
minor axis of the ellipse before a approaches . 


p= fw dy’ = Ff Nea ay SE 
Te rea ie en ey fae ba ieee apie: 


Asa = o, let P = o in such a way that P/a remains equal to a 
finite constant. Then the value of P from (68) will in the limit be the 
force per unit length of the strip of pressure. Then fora = « 


By Bs ae 2P yf? 
ape es Hota ona 
ex 27ab \: b> orb : Z (69) 


Since (Appendix B) at a distance r the potential of an infinitely long 
straight line of constant linear density P’ is — P’ log r? 


b 
o= aa P’ log r2dy’ (70) 
Sh. 


where r= (y — y’)? 4+ 22 (71) 
and where P’ is given by (69). 


By differentiation under the integral sign 


d¢ 0d 0° x * Od 


d¢ b pr , ao b pl b pr 
eee ie) ee -2 f aay +40 fi a dy’ (73) 


Now for points on the z axis, i.e., for y = 0, we have by (69) 


pe oP ey Z 
~ ay = =f 3 ee hop 
-b Tv aie b2 
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and letting y’ = b sin 6 
ie : 2P 3 b2 cos?6 dé 
Let ee ae, -3 b?sin?9+4+ 2? 


2 b2 2 b2 2 a 
ae | - @ + Nee are tan eure tno) |? 
2 


arb? 


i 2P (Vb? + 22-2) 


= (74) 
wae ee yi as 
ie —b ip 7 ab =i} fe a b? Y 
bu 
and letting y’ = Vit wv 
1+ uv? 
pt 2P ‘ du 
dy’ = a [ (b? + 22)u? + 2? | 2 


» 7 


a : oe A 
ie bate + 22)u? + Z| 3 223 Vo: ies are tan ; 


00 


eee 
= 2e\b? + 2? ey 
Substituting the integral evaluations from (74) and (75) in (738) 
og *( eu ) 
Gy a Eek a 6 
= wl bese =e (76) 
a2 AP pe Aue 
a wel Vo* ae :) MF Vo? ae ae 


By substituting values of the derivatives of ¢ from (72), (76) and (77) 
into (9) 


4oP —— 
oe eee (vir cera :) (78) 
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2P 
Lipa ie toe 79 
f r\b + 2? ED 


To get Y, it is noted that the stress at every point is a function only 
of y and z and hence the strain is also a function only of y and z. 


Thus du Aue ov Ow 
go Ra Oe 


Hence from (8) we have: 


dv dw 
Y,+ Z, = 2dA + 2y By ee = 2(A + w)A (80) 
Xe =A 2 a 


and hence by means of (4) 


Xe = ¢(%, + Z, 
or Xz 
Yy oe oe = 45 (82) 


Substituting in (82) the value of X, from (78) and that of Z, from (79) 
oP (Vb? + 22 — 2)? 


ye ai ae ey 
By means of (54) and (68) 
oe 2rb2 a idle 8 ab? 
WGaekP © oe nae or ee 


Hence with the values of A, X., Y,, Z. from (84), (78), (83), and (79) 
respectively, we have 


RX. eee ed 

; 20 (V1 + (2/b)? — z/b) 
nee (Vi + @/b)? — 2/b)? 
b NL + (/b)2 
AZ ee 
b Vi + G/b)? 


(85) 
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The stress components are thus algebraic functions of z/b and it is 
readily found by the usual process that |Z, — ¥;| will be a maximum 


5-1 
for 2/b = yee = 0.7861. For z/b = 0.7861 we have from 
(85) : 


OG 
= — 0.97186 
b 
MY 
5 = — 0.1856 (86) 
AZ, 
= — 0.7861 
b 
A 
and b (Z, — Yy) = —0.6005 (87) 


Ill. Merruop or Souution 


13. Notation and Equations for Definite Problems.—For the con- 
venience of those readers who find the mathematical derivation just 
presented to be difficult reading, below are given the essential nota- 
tion and equations to be used in the solution of definite problems. 
The subscripts 1 and 2 are used to indicate one or the other of the 
two bodies in contact. The following notation is used: 

E,, E, Moduli of elasticity for the two bodies in contact. 
o1, 2 Poisson’s ratios for the two bodies. 
B Mean of reciprocals of the radii of curvature of the two 


1 i] 
bodies in the y direction. B= 4% (= + =) 
tie © Ate 
A Mean of reciprocals of the radii of curvature of the two 
1 I 
bodies in the z direction. A = 4% (= + 5). Note that 
Wiha ed es 
B/A is always equal to or greater than unity; this fact 


determines the x and y directions. The ratio B/A is 
related to the value of k as shown by Equation (52). 


1 — o;? 1 — =| ] § 
A= | a + E, TORE (See Equation (55)) 
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Note that this quantity A is a term which takes account of the 
elastic properties of the two bodies in contact and also of the shape of 
the two bodies at the point of contact. If the two bodies in contact 
are similar, such that 01 = o. = cand H; = Ey = EH, then the above 
reduces to 

2(1 — o?) 
A = RGB) 


a is the semi-major axis of the ellipse of contact. 
b is the semi-minor axis of the ellipse of contact. 
k is the ratio b/a. 

ki = 1 — ke? 

P is the total pressure between the two bodies. 


{ 3 APE (k’) 
a= 
Qrk? 

where E(k’) is a definite elliptic integral. (See Equation (41)). 

X,, Y,, Z, Normal components of stress across planes perpendicular 
to the x, y, and z directions through the point under con- 
sideration. Attention is called to the fact that the stress 
components considered in this bulletin are for points on 
the z axis at various relative depths, z/a, from the contact 
surface (see page 17). 


A A A 
The values of * ee - Y,, and a Z, are given by Equation (59) 


which will not be repeated here. The methods of solving those equa- 
tions are discussed below. 


14. Stress-Component Factors.—As suggested on a previous page, 
the substitution of numerical values for the purpose of solving 
Equations (59) for the stress components is practicable only by the 
choice of certain values of k which will permit obtaining the values of 
the elliptic integrals involved from published tables. The values of k 
which can be substituted in this way probably will not be in agree- 
ment with the numerical values for which a solution is desired, since 
for convenience in the substitution it is desirable to take k = cos 6, 
where 6 may vary from zero to 90 deg. If Hancock’s table of Elliptic 
Integrals is used, 6 should be some multiple of 5 deg. It will be found 
much easier to substitute in these equations if ¢ is taken in degrees 
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2 
such that the values of aA = cot ¢) will be approximately those de- 


sired, instead of assuming values of z/a and solving for yg. In order to 
facilitate the solution of practical problems, the right hand members 
of equations (59) have been evaluated for values of k ranging from 
1 to 0.2588 and for z/a ranging from 0 to 4. (In some eases the values 
of z/a are from zero to one). Poisson’s ratio, o, is assumed as 14. 
The computed values are given in Tables 1 and 2. In order to obtain 
the stress components, the factors in these tables must be multiplied 


a 
by a? determined from Equations (54) and (57). In Fig. 1 have 


been plotted for a value of k of 0.86603 the values of the factors for 
the three stress components from Table 1, using the values of z/a as 
abscissas. To avoid confusion the individual points are not indicated. 
In this figure is also shown the curve for the maximum difference of 
stress components, (Z, — Y,). It is to be noted that the values of 
the individual stress components are a maximum at the surface 
(z/a = 0) and decrease with increase in distance from the surface of 
contact. However, the values of (Z, — Y,) increase from a small 
value at the surface to a maximum at a point below the surface where 
z/a = 0.44 (for k = 0.86603). So far as tendency to cause damage 
is concerned, it is evident that for k = 0.86603 the greatest value of 
the difference of stress components, called (,Z, — iY,), is equal to 


a 
— 0.384 A and that it occurs at a point below the surface z:/a = 0.44, 


as marked on the diagram. These curves are shown to illustrate the 
method used in obtaining curves for the maximum differences of 
stress components in Fig. 2, and to indicate how the stress compon- 
ents vary with the distance from the contact surface. 

In order to determine the greatest value of the difference of stress 
components (1Z, — 1Y,) for the other values of k, curves similar to 
those of Fig. 1 were plotted from quantities given in Tables 1 and 2 
and the values thus determined for each value of k are given in the 


A A 
lower part of Table 3. The corresponding values of 7 1X5, - eae 


A 21 


— ,Z,, as well as of a? ure also included. The stress component 
a 


a . 
factors ( Xz, etc.) at the surface of contact, obtained from Tables 1 
i 
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TABLE 3 
SuMMARY OF Factors FOR CALCULATIONS 


These factors are to be used for calculating stress-components, size of contact areas, and depth to 
point of maximum stress difference. Poisson’s ratio = Ye 


Values of k 

1 0.86603 | 0.70711 | 0.57358 | 0.42262 | 0.34202 | 0.25882 
BY Apes ete 1.0000 | 1.2409| 1.6829] 2.3083] 3.6746| 5.0831 | 7.8625 
eT eee 1.3333 | 1.0704 | 0.7753 | 0.5474 | 0.3214 | 0.2191 | 0.1303 
a 
Ts Coes —0.4775 | —0.4320 | —0.3702 | —0.3109 | —0.2355 | —0.1919 | —0.1451 
ah eee ~0.4775 | —0.4532 | —0.4151 | —0.3726 | —0.3092 | —0.2669 | —0.2160 
2 oh eee ee ~0.6366 | —0.5901 | —0.5235 | —0.4557 | —0.3631 | —0.3058 | —0.2407 
“ (oz, — oXz}..) —0.1591 | —0.1581 | 0.1533 | -0.1448 | —0.1276 | —0.1139 | —0.0956 
» {o2, — o¥y}--| —0.1501 | —0.1369 | —0.1084 | —0.0830 | —0.0539 | —0.0390 | —0.0247 
* 7 aa ee —0.1130 | —0.105 | —0.100 | —0.090 | —o0.082 | —0.070 | —0.056 
Patras tant —0.1130 | —0.097 | —0.084 | —o0.064 | —0.052 | —0.041 | —o0.028 
* Ay AAR —0.5230 | —0.481 | —0.427 | —0.367 | —0.293 | —0.243 | —0.187 
*. {:22—-1¥y}..| -0.410 | —0.384 | —0.343 | -0.303 | —0.241 | -0.202 | —0.159 
Freee he de 0.467 0.440 | 0.390 0.350 | 0.280 0.240 0.200 


and 2 for z/a = 0, have been included in this table, as have also the 


AP 
computed values of ae (from Equation (54)) and of B/A (from 
Equation (52)). 


15. Method of Substitution.—Using B/A as abscissas, the values 
INP 21 A A 
GZ.— 1Yy) and 7 (oZ2 — oXz) from Table 3 have been 


plotted in Fig. 2. This diagram (or a similar diagram plotted to a 
larger scale) permits the determination of these factors for values 
B/A other than those given in the table. In making use of this dia- 
gram for the solution of a numerical problem, the first quantity to be 
determined from the conditions of the problem is the value of B/A, 


1 1 
where A = OR,” and B= OR, (for crossed cylinders having radii 


of Ri and Re, respectively, and with axes perpendicular). Then for 


‘ AP. 
this B/A, the value of Ve is taken from the diagram, which permits 
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solving for the value of a for given values of P and A. Having cal- 
culated a and A, (.Z, — 1Y,) may be computed by multiplying the 
A a 

factor for - (iZ, — 1Y,), taken from the diagram, by A: Within the 
limits covered in this diagram, Z, — Y,is negative in sign. The depth 
(in inches) to the point of maximum stress difference may be obtained 
from the curve for z:/a, the value of a already having been deter- 
mined. The semi-minor axis of the ellipse of pressure b is equal to ka 
and if desired the pressure area may be calculated as rab. 


16. Surface Stresses.—If the stress components at the surface of 
contact are desired it will be necessary to construct a diagram similar 
to that shown in Fig. 2, using for the purpose the values of the surface 
stress-component factors in the upper part of Table 3. While the 
stress components at the surface will be higher than those in the 
interior, the stress differences will be less. To show the comparison, 
a curve for the maximum difference of stress components at the sur- 


A 
face, as (oZz. — oXz), 18 included in Fig. 2. Attention is called to the 


fact that at the surface the maximum stress difference is (Z, — Xz), 
while in the interior it is (Z, — Yy,). 


17. Numerical Example.—As an illustration of the method to be 
used in substituting definite values in the mathematical solution, a 
numerical problem will be solved. The radii of curvature and shear- 
ing stress used in this problem are those used in the check of the 
mathematical analysis, as described in detail later. Assume crossed 
cylinders of steel having axes perpendicular to each other, the radius 
of one being 15.8 in. and of the other 12 in., Poisson’s ratio ¢ = 4, 
and the modulus of elasticity H = 30 000 000 lb. per sq. in. 

For the problem under consideration 


1 1 
A aN 15 o UU: 
1 1 
BS QR, ~ 9x 12 = 0.0417 
af 0.0417 Cae 
B/A = On0346e. <5 


B+A = 0.0733 
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In Fig. 2 for B/A = 1.320, from the corresponding curves the fol- 
lowing factors are found: 


k = 0.838 
AP 
— = 1.010 
a 
al 
— = 0.428 
a 
i ’ 
ih GZ, a: ig) = —(0.375 


From Equation (55), assuming o1 = o2 = o; H; = Ey = E, we have: 


pees) 2(1 — 0.0625) 
A = (A+B) ~ 30000 000 X 0.0733 


(ee, 0.00040 YP 
MEA [SHRINE wage 


a 3 
{= 11090 yP 


The maximum difference of stress components is 


= 8.52107" 


a 2 a 
ihe —1iYy = —0.375 rcs — 4160 y? 


In order to calculate the load required to produce a shearing stress of 
19 200 lb. per sq. in., corresponding to a stress difference of 38 400 lb. 


3 
per sq. in., we have 38 400 = 4160 \P 


from which 


(= 400 


3 
aa) = ‘7AXD Moy. 


For calculating the area of contact, 


See 3 
a = 0.00946 \? = 0.00946 4786 = 0.0873 in. 


b = ka = 0.838 X 0.0873 = 0.0731 in. 
Area of contact = rab = 0.0200 sq. in. 
Depth to point of maximum stress difference, 

Z1 = 0.4284 = 0.428 X 0.0873 = 0.087 in. 
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786 
= 39 300 lb. per sq. in. 


verage surtace pressure 0.0200 


The maximum surface stress must be calculated from the values 


mA 
of "7 Zo in Table 3. 


Additional numerical examples will be found in Appendix C. 


18. Special Case.—It may be of value to indicate one further point 
concerning the method to be used when one (or both) of the bodies 
has more than one radius of curvature at the point of contact. 

It will be remembered that in the definition of the quantities 
A and B, it is assumed that the shape of the bodies in the neighbor- 
hood of the surface of contact may, by a proper choice of axes, be 
represented by a quadratic of the form 

2 = Ax? + By? (See page 14) 
For crossed cylinders with axes perpendicular, as already indicated, 
il 1 
CES TDR BS 2R2 
where R; and R» represent the radii of curvature of the two cylinders 
in the x and y directions, respectively. But if one of the bodies has 
one radius of curvature in the a direction, Ri, and another radius in 
the y direction, R’1, while the second body has only, say, the radius 
R, in the z direction (in other words, R’, = «) then 


1 1 
A=4(Etz) and 


- 1 1 1 
B=% R’; = Pe a [Pore QR, 
That is, A is the mean of the reciprocals of the radii of curvature in 
the x direction while B is the same function for the y direction. After 
A and B have been calculated, the remainder of the computations 
may be made as just described for crossed cylinders. It should be 
remembered that the radius of curvature is to be considered as nega- 
tive if the center of curvature lies outside of the body. 


f 
IV. CuHeck or MATHEMATICAL SOLUTION 


19. Need for Checking. —In applying the mathematical calcula- 
tions of stress components to certain practical problems, the existence 
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of such high shearing stresses was indicated that it was thought desir- 
able to attempt to obtain a physical check of the correctness of the 
mathematical assumptions and computations. 

In the paper previously referred to, Hertz described the method 
used by him for verifying his analysis. For this purpose he deter- 
mined the areas of contact between a spherical glass lens and a glass 
plate at various loads, and also between crossed cylindrical glass rods. 
Since the measured areas showed good agreement with those calcu- 
lated, he took this as a check also of the stress components. 

Because of the difficulty of measuring accurately the areas of con- 
tact, much more difficult for metals than for glass, and also in order 
to obtain a more definite check of the calculated maximum shearing 
stresses, it was decided to make use of the Fry strain-etch method* 
for determining when a material has been stressed beyond the yield 
point. 


20. The Fry Strain-Etch Method.—By the Fry method of etching 
(which is applicable only to certain low carbon steels), after the 
material has been subjected to a stress exceeding the yield point, by 
bending, direct loading, or other means, it is heated for half an hour at 
a temperature of approximately 400 deg. F. After heating, the piece 
may be cut at any section which it is desired to examine and on etch- 
ing with a solution of cupric chloride in hydrochloric acid and water 
it will be found that those regions which have been overstressed will 
etch darker than the remainder of the section. According to Fry’s 
method, etching was accomplished by immersing the specimens in the 
etching solution for as long as seemed necessary, followed by rubbing 
with copper chloride crystals wet with the solution; but it was found 
that etching could be accomplished very readily and more rapidly by 
rubbing the surface with a swab made by twisting cotton around a 
stick, the swab being dipped in the etching solution and in copper 
chloride crystals from time to time. After a satisfactory etch has been 
obtained the surface is cleaned with alcohol containing a small amount 
of hydrochloric acid. The acid is used to prevent copper from being 
deposited on the surface. 


21. Special Materials and Equipment.—For this method of check- 
ing the mathematical analysis there was obtained a bar of mild steel 
24 in. square and 6 feet long, which by preliminary tests it was found 


*“Strain Figures in Mild Steel Evidenced by a New Etching Method,” Ad. Fry, Stahl und Eisen 
Vol. 41, No. 32, Aug. 11, 1921, pp. 1093-97. A condensed translation of this paper may be found in the 
Tron Age, Dec. 1, 1921, p. 1401. Additional work along these lines is given in a paper by J. Dudley 


‘ ee eee Detection in Mild Steel by Special Etching,’’ Journal of the Iron and Steel Institute, 
ol. 5 
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could be used. It is a peculiar fact that all mild steel will not give 
strain-etch figures, but the reason for this is not known. Therefore, 
in order to obtain a material which would strain-etch, six bars were 
obtained from different sources and it was found that one of them was 
satisfactory. The average physical properties of this material were: 


mensilerstrength. .. 54.64.24 +s 47 400 lb. per sq. in. 
ensiletyield pointe. o.y.<.. + 32/00 1b. per sae in 
Compressive yield point........ 28 800 lb. per sq. in. 
“orsional yield point. .2... 19 200 lb. per sq. in. 


(solid specimen) 
So far as these tests are concerned only the yield points are of interest. 
Two types of loading were used in checking the mathematics. The 
first consisted in loading a plane specimen block through a cylindrical 
loading block; the second, in loading a cylindrical specimen block 
through the same loading block, the axes of the cylinders being at 
right angles to each other. 

The loading block consisted of a segment of a cylinder whose 
radius was 15.8 in. The length of the block in the direction of the 
axis of the cylinder was 214 in., its radial thickness 334 in., and its 
width 6 in. The block was machined to size, and then the cylindrical 
surface was machined, the radius used being the largest which could 
be obtained with the available equipment. After case-carburizing 
and quenching, to prevent over-stress of this block during loading, 
the cylindrical surface was ground to as accurate a finish as possible. 

The specimen blocks for the tests of stresses between a cylinder 
and a plane were machined from the steel already described. These 
blocks were made 214 in. square and 12 in. long, one face being fin- 
ished as accurately as possible to a true plane surface. For checking 
the accuracy of the surface, the blocks were coated with prussian blue 
and when tested with a straight edge showed practically perfect 
contact. For the crossed-cylinder tests, a cross-radius of 12 in. was 
machined on one of the square blocks just described. Considerable 
difficulty was encountered in obtaining the degree of accuracy of 
finish desired, both for the loading block and the specimen blocks, 
but it is felt that the surfaces finally obtained were satisfactory. 


22. Tests of Cylinder on a Plane.—In the case of a cylinder on a 
plane, the calculation of stress components for a given load and a 
given radius of a cylinder is relatively simple. By referring to 


2PA 
Equation (84), page 28, it is found that b? = = where b is one-half 
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of the width of contact, P is the load in pounds per inch of length of 
2(1— 0%) 4R(1 — o?) 


the cylinder, and A = 1 = E . For steel, assuming 
POR 
Poisson’s ratio, c, as 0.25 and E as 30 000 000 Ib. per sq. in., we have 
4R(1 — 1/16) 
~ 30.000 000 
R 
me My 
2 Pi a 
= V— —___ = 0.000282 V PR 
Then b Bee 0.000282 Y 


— 0282 SONS 10° 
n 0.00 x OS 


2256 \ a 
R 


The maximum difference of stress components, from Equation (87), 


I 


I 


b 
is, eal, = — 0.6005 > 


— 1355 y! 
R 


For the loading block used in these tests, R = 15.8 in., and hence 


v 


ll 


(Zee a ea IAN 


Remembering that the shearing stress, S, is one-half of the stress 
difference, we have 


341 j= = 
S = (2, — 1¥,)/2 = > VP = 170.5 VP 
The load required, then, to produce an intensity of shearing stress 


equal to the shearing yield point of the steel used will be 


= (= 200 


2 
= = 12 700 lb. per in. of length 
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LQ 9 
Ss 9 

Ss Q 8 
Nv 2 . 


Fig. 3. Srrain-EtcHep SpecIMEN Biock 


For this load, the width of contact is 


2b 


2 X 0.000282 V/12 700 X 15.8 
0.000564 X 448 
= ().252 in. 
Depth to point of maximum shearing stress 
2, = 0.786b = 0.786 X 0.126 = 0.099 in. 

In order to check the load calculated as producing a shearing yield 
point stress, the specimen block prepared from the strain-etch steel 
was loaded at different points along its length with five loads as 
follows: 


Load 
Load No. lb. per in. Remarks 
i Ls 8 3 oh eaten Eas aOR Ree 12 000 No strain lines 
DAE sts COE A ST TE RRR INNS 12 450 No strain lines 
Ne A Berk OE UE ORES OR CTE DRIAL OC MSE 12 930 Strain lines 
Cio des deege gO mana Os ay RAN Don hai, Net reir 15 100 Strain lines 
1 ae OR Be eee sn st eee re eT 17 800 Strain lines 


It will be noted that two of these loads are lower than the 12 700 lb. 
per in. calculated above, while three are higher. The specimen block, 
after heating as already described, was split longitudinally, the sec- 
tion finished smooth, and then etched. It was found that for the load 
of 12 930 lb. per in. and for all higher loads the sections showed 
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strain-etch patterns, while for the two loads lower than 12 700 lb. 
per in. the sections showed no evidence of overstress. 

Figure 3 is from a photograph of sections under loads 1, 4, and 5 
and shows the kind of strain-etch patterns obtained. The magnitudes 
and points of application of the loads are shown in the figure. It 
should be noted that the strain lines practically disappear at the sur- 
face and that the lines are at an angle of approximately 45 deg. at the 
distance below the surface (0.1 in.) where the maximum shearing 
stress is expected. Only small portions of the specimen block are 
shown in this figure, since it was found desirable to cut the specimen 
block into relatively small pieces for convenience in etching. 

This same procedure was repeated with a second specimen block 
and again the strain-etch figures were not obtained until the load was 
such that the calculated shearing yield point stress of the material 
was exceeded. The accuracy of the determination from the loads 
appears to be on the order of plus or minus two per cent. This is 
probably within the limits of accuracy of the determination of the 
shearing yield point. 

It may be worth while to compute the maximum intensity of the 
direct compressive stress (o9Z,) at the surface under the load of 
12 700 lb. per in. . 

From Equation (85) we have 


b 
Zz SS A 
We have already found that 
uy = 2256 \ a 
A R 
Hence: 
oZ2 = —2256 qe —63 900 lb. per sq. in. 


In general, for this type of loading (cylinder on a plane), the maxi- 


mum stress at the surface is equal to s times the average stress or 
AP ‘ 

a 2b° 

While the compressive yield point of this material was only 


28 800 lb. per sq. in., yielding in direct compression did not occur at 
the surface because of the high stress components at right angles to 


0Zz = 
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this. From Equations (85) it may be determined that at the surface 
oY, is equal to oZ., and oX,, is one-half as great, if Poisson’s ratio is 
taken as 44. That is, the maximum shearing stress at the surface is 


63 900 
only 470 approximately 16 000 lb. per sq. in., as compared to 


19 200 lb. per sq. in. at a point in the interior. 


23. Tests of Crossed Cylinders.—In checking the mathematical 
analysis for stresses within the area of contact of crossed cylinders, 
the same procedure was followed as has just been described for a 
cylinder on a plane. After calculating the load required to produce a 
shearing stress equal to the shearing yield point of the material a 
series of loads extending below and above the calculated load was 
applied, and after heating, the block was split, etched, and examined 
for evidence of overstress. 

In this case the mathematical computations are much more diffi- 
cult than for the loading previously considered, unless the empirical 
equation discussed in the following pages is used. Substituting the 
conditions of this loading in that equation 


22 00m ee 
1Zz coal Tey Se See 7 ear aee 2% 
Cefityy joa R;' 
(4S 0.271 3 
or P = |2~x19200\ 12 ) A) (12)? 
23 500 


= 9.233 = 786 lb. 


It is seen that this value of the load agrees with the value previously 
obtained for these same conditions, using the exact method with 


A 
graphical interpolation for a GZ. — 1Y,). 


After loading the specimen block with a series of loads ranging 
from below the calculated load to well above it, the block was split 
through the centers of the loaded areas and on etching was examined 
for evidence of overstress. It was found that for this condition of 
loading the strain-etch lines were not clear and distinct as for the 
loading first used. Although the lines were visible during the etching 
on those sections where the load was such that the computed shearing 
yield point was exceeded, it was found that on cleaning the surface to 
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prevent rusting or the deposition of copper, the lines became so 
indistinct that they could not be photographed. Much effort was. 
spent in the attempt to produce etch-lines which could be photo- 
graphed, but without success. However, it may be stated definitely 
that the mathematical computation for crossed cylinders has been 
verified within reasonable limits by means of the strain-etch method. 


24. Approximate Equation—Crossed Cylinders.—As has been indi- 
cated, the direct use of the mathematical analysis for the solution of a 
definite problem is a long and tedious process. The most satisfactory 
general method of substituting definite numerical values is by the use 
of a diagram such as Fig. 2. However, for the case of crossed cylinders. 
it is possible to derive an approximate empirical equation which per- 
mits an easy solution for values of the maximum difference of stress. 
components. Using the two equations ; 


AP IAP 
SP Tyme a 
a D 


where D is a constant (see second line, Table 3) for a given value of 
B/A, and 
2(1 — o?) 
~ E(A + B) 
assuming that both bodies are steel and that o = 144, E = 30 000 000 


lb. per sq. in., and remembering that R; = (B/A)Rz, it can readily 
be determined that 


a= 40 000 Pe 
A | B/A i R,” 
B/A + 


Referring to Table 3 it will be noted that for a given B/A, the 
value of (,Z, — iYy) is equal to a numerical factor multiplied by =. 


If then for each value of B/A in Table 3 we multiply the right-hand 
side of the preceding equation by this factor, and substitute the cor- 
responding numerical values of B/A and D, it will be found that the 
maximum stress differences may be represented by the equation 
Pp 
12, = mee =C€ R® 
where C is some function whose form is to be determined, but where 
C has a definite numerical value for a given B/A. If the values of C 
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obtained as just described for values of B/A in Table 3 are plotted on 
log paper using the corresponding values of B/A as ordinates, it will 
be found that the points lie on a straight line giving the equation 


oe 23 500 
in? (B/A)?271 
Hence we may write: 
23500 P* 
14, — eg = (R1/ Ry) 27 R.* (a) 


For convenience in solving with a slide rule, this may be changed to a 
slightly less accurate approximation: 


Wea 23 000 P* . 
V2 ea (Ry /Re) R,” ( ) 


: R 
Equation (b) gives values about two per cent low for — = i 
i 2 


R 
and about 21% per cent high for a = 7.86. The error is less for 
2 


B 
intermediate values. So long as of is kept within the limits found in 


Table 3 these equations may be used in solving for maximum differ- 
ence of stress components. While the equations appear complicated 
they may be solved with a slide rule, although (b) is the more easily 
handled. In Fig. 4, Equation (a) has been plotted in the form of a 
nomographic chart which may be used for solving for any one term 
if the other three are given. In using it the intersection with the 
“support” (ungraduated line) of a line joining known values on the 


R : aay : 
lines for P and — must also be on the line joining corresponding 
2 


points on the ,Z, — 1Y, and R, line. As an example, if P = 10 000 
R R ; 
Blb:, — = 1, R, = 12 in., first connect =a with P = 10 000, then 
Ro Re 
from Ry = 12 draw a line through the point where the first line cuts 
the “support.” This line, extended, intersects the 17, — iY, line in 
the value desired, 96 500 lb. per sq. in. 
As in the case of the empirical mathematical equation, the use of 


R ie ake ; : 
this chart for values of es outside of the limits indicated is not advis- 
2 


able, for the errors involved may be considerable. 
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Values of Ss a, and k, may be more readily obtained from the 
a 


corresponding curves of Fig. 2 than by the use of any empirical equa- 
tions which might be derived. 


V. Errect or Various Factors oN CALCULATED STRESSES 


25. Significance of Calculated Stresses.—By the use of the method 
of calculation presented in this bulletin, for many conditions it is 
possible to compute the stress components and maximum shearing 
stresses produced when one elastic body is pressed against another. 
In general, the conditions of loading will probably be such that during 
the useful life of the member it will be subjected to a large number of 
applications of the load. Such loading produces a shearing stress 
which varies from zero to a maximum for each application, and hence 
it is evident that for materials used under this kind of loading the 
fatigue endurance limits of the materials when subjected to a shearing 
(torsion) stress varying from zero to a maximum is one of great im- 
portance. There is not a great amount of information available 
concerning the endurance properties of materials subjected to this 
kind of stress, and hence the interpretation of the significance or 
seriousness of a particular calculated shearing stress is open to doubt. 
An important question which cannot be settled at this time is 
whether or not a repeated shearing stress having a maximum value 
in the interior of a body will cause a fatigue crack to form and spread 
under this condition as readily as in the case of a torsion fatigue test 
specimen where the maximum stress is at the surface. Some prelim- 
inary tests made in the Fatigue of Metals Laboratory have a bearing 
on this question, and the results seem to indicate that cracking and 
rupture may not start in the interior of a body when the computed 
shearing stress is somewhat higher than the torsion endurance limit 
determined in the ordinary way. This needs considerable further 
study, especially since it is difficult to make tests of this nature with-_ 
out having a change occur in the curvature of the specimen due to 
wear, and hence a quite appreciable change in calculated stresses. 

Another question arising in this connection is the possibility of 
work hardening of the material due to overstress, thus causing a 
change in its physical properties. The residual internal stresses due 
to cold work and rolling action may also have an effect both on the 
amount of shearing stress produced at the point of contact of the two 
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Fic. 4. NomocrapHic CHArtT For APPROXIMATE EQUATION 


bodies and possibly also on the endurance properties of the material. 
In fact, it would seem that under some circumstances residual internal 
compressive stresses due to cold work might so affect the shearing 
stress as to cause the actual maximum shearing stress to be at right 
angles to the direction which would be expected without taking ac- 
count of the stress due to cold work. That is, the maximum stress 


46 ILLINOIS ENGINEERING BXPERIMENT STATION 


difference might change to Z, — X, instead of Z. — Yy,. This would 
depend to some extent on the relative radii of curvature of the two 
bodies, and for this to occur the internal stress due to cold work 
would have to be greater than X, — Yy. 


26. Crossed Cylinders.—In order to study the effect of the various 
factors entering into the equation for stress differences the empirical 
equation 

23500  P* 
(Ri/R2) 0.271 R.* 


12, a hs = 


can be changed to the form 


23.500 P% 


R,°?27 Jip ORte 


iZz aa ay = 


From this it is seen that the maximum stress difference varies as the 
cube root of the load, that it varies inversely as R,°-?7! and inversely 
as R,°-*%, In this case Ri > Re, and hence it is seen that the effect 
of changing the smaller radius a given amount is greater than that 


for a corresponding change in Ry. Again, attention is called to the. 


R 
fact that these conclusions are valid only for values of — between 


Re 
the limits of 1 and 8. 


27. Cylinder on a Plane.—For the condition of loading of a cyl- 
inder on a plane, the maximum stress difference varies directly with 
the square root of the load per inch of length of cylinder and inversely 
as the square root of the radius of curvature of the cylinder. That is, 
the stress does not vary so long as the ratio of the load to the radius 
of curvature is a constant. 


28. Effect of Poisson’s Ratio.—In all calculations in this bulletin, 
Poisson’s ratio has been assumed as 4. ‘The effect of changing 
Poisson’s ratio to some other value is to change the stress component 

A A 
factors 5 XxX, and - Y,in Tables 1 and 2. In order to determine the 


A 
effect on if GZ, — iYy) of assuming this ratio as 4 instead of 14, 


the values of the stress-difference factors have been calculated for 
o = 14 for the various values of k. It is found that the values of 
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A 
m (:Z, — iYy) are reduced, when o is changed from 14 to Y%, in the 


ratio shown in the second column below. 


Stress-difference 


_ Value of k Raves i Net Effect 
a= hy 
ILS OOO Se Feresscieueesnoe cere eae 0.94 0.97 
ORSGG Mee rs epee anole eat 3 0.95 0.98 
ORLO Teor Cinta tee eas 0.95 0.98 
ORS Sues etches toes 0.95 0.98 
(0) Pe a See ie ea neg ete 0.96 0.99 
RGA Mites tics rota em oan tees 0.96 0.99 
OR OSaer oteeaa heir den 0.96 0.99 


a 
However, we find that for crossed cylinders the quantity 7 is increased 


by 3.5 per cent when oa is changed from 14 to 4. The net effect then 
is to change the ratios to the values given in the last column. For 
other than crossed cylinders it will be necessary to compute the effect 


a 
on 5 of a change in Poisson’s ratio. If desired, calculated stress 


difference may be corrected by the factors just given, but it should 
be noted that in using the tables or diagrams for making calculations 


a 
the value of K must be based on o = 14. The broken line in Fig. 1 


shows the stress differences for k = 0.86603 and o = 4. 
At the surface the change in Poisson’s ratio has a greater effect. 


A A 
At this point the values of = (oZ. — oY y) and of z (oZ2 — oXz) for 


o = 4 are two-thirds as large as for o = 4, for all values of k. 

For the case of a cylinder on a plane, there is no difficulty in direct 
substitution of any value of Poisson’s ratio desired, since the equations 
are readily solved. 


VI. SumMMarRy AND CONCLUSIONS 


29. Conclusions.—As the result of the investigation the following 
conclusions were drawn: 

(1) By using the methods of solution presented in this bulletin, 
it is possible to compute the principal stress components within the 
area of contact, both at the surface and at points below the surface, 
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when one elastic body is pressed against another. It is also possible 
to compute the maximum shearing stress, due to the difference of 
stress components, and to determine the distance from the contact 
surface to the point of maximum shear, as well as the area of con- 
tact. 

(2) The accuracy of the calculations of maximum shearing stresses 
‘in the interior of a body has been verified by a strain-etch method. 

(3) For the case of crossed steel cylinders, if Poisson’s ratio = \4, 
the maximum stress difference is given by the following empirical 
equation: 


J 


23 500 ee 
(Ry / Re)?" R,# 
This maximum stress difference (which is equal to twice the shearing 
stress) occurs at a point below the surface of contact. This equation 
may only be used for values of Ri/R2 between the limits of 1 and 8, 
since calculations have not been carried beyond those limits. R, is 
the radius of curvature (expressed in inches) of the larger cylinder, 
and RF, is the radius of curvature of the smaller. P is the total pressure 
between the two bodies, in pounds. 


12, ion we = 


The exact solution for this condition and also for more general. 


conditions is found in the text. 
(4) For the case of a steel cylinder on a plane, the maximum 


ihe 
shearing stress equals 677 @=, the width of contact is 2b = 


0.000564/PR , and the depth to point of maximum shearing stress is 
2: = 0.786b = 0. 000222\/PR . The maximum intensity of surface 


Lead os 
pressure is equal to So 2256 Rp That is, the maximum 


surface pressure is 27 per cent higher than the average. In this case 
P is the load per inch of length of the cylinder whose radius is R. 

(5) The mathematical analysis is based on the assumption that 
the bodies in contact are homogeneous and isotropic. While it is 
well known that shapes and sections rolled from metals do not con- 
form exactly to this assumption, it is probable that the errors in- 
volved will not be large. The question as to the seriousness of a 
given calculated shearing stress is one which can not be answered 
definitely until more information is obtained on the endurance 
limits of materials as commercially produced, when subjected to 
repeated shear, and also the effect of cold work and the residual 
stresses caused by cold work on the physical properties and behavior 
of such materials. 
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APPENDIX A 


REDUCTION OF A QuapRATIC ForM To Ax? + By? 


The expression A’é,? + 2H’&.y, + B’n,? must be reduced to Ax? 
+ By? by a rotation of axes through angle 6. This rotation is 
effected by making the substitutions 


HS is ea (1) 
m1 = «sind + y cosé 
Thus we have after collecting terms 


A’t? + 2H’ €in: + B’n? = 22(A’cos? 6 + 2H’sin 8 cos 6 
DUB icin: 0) 


+2zy[(B’ — A’)sin6é cosé@ + H’'(cos?6 — sin? 6)] (2) 
+y?(A’sin?@ — 2H’ sin @ cosé + B’cos? 6) 
Equating this to Ax? + By?, it is found that 
A = A’cos?6 + 2H’sin 6 cosé + B’sin?6 
0 = (B’ — A’)sinécos6 + H'(cos?6 — sin? 6) (3) 
B = A'sin?@ — 2H’sin 6 cos6é + B’cos?6 
It may be verified that 
A Bis AB" (4) 
AB = A’'B’ — H’? + [(B’ — A’)sin@ cos@ 
+ H'(cos?6 — sin?6@)]? = A’B’ — H” (5) 


Now the quadratic equation x2 — (A + B)x + AB = 0 has roots 
A, B. With the values of the coefficients just found, the quadratic 
equation with roots A, B becomes 


a? — (A’ + B’)rz + A’B’ — H”? =0 (6) 


From (12) we have (see Chapter II): 
2(A’ =e B’) = Pil a= p'1 te pe a p's (7) 


4A'B’ = [pi + p2 — (p2 — p’s)sin? el[p’1 + p's + (p2 — p’2) sin? ¢] 
= (p1+ po)(p'1 + p’2) + (02 — p’2)(p1 — p'1 + p2 — p'2) Sin? € 
— (po — p’s)*sin‘ e 
4H’? = (po — pi)? sin? e cos? e = (pz — p’2)?sin? € — (p2 — p’2)’sin* € 
4(A’B’ — H') = (pi + p2)(p'1 + p’2) + (01 — p'1)(p2 — #’2)sin?e (8) 
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Substituting the values of A’ + B’ from (7) and A’B’ — H”? from 
(8) into (6) and multiplying the resulting equation by 4, we get 
427) — 22(p1 -F Deka (5 ar p’2) os (pi = p2)(p'1 ch p’2) 
+ (pi — p’1)(p2 — p’2)sin? e = 0 
or (2% — pi: — p2)(2e — p': — p’2) + (pi — p's) (p2 — p’2) sin? € = 0 (9) 
The two roots of (9) are the values of A and B desired. 


APPENDIX B 


POTENTIAL OF AN INFINITELY LONG STRAIGHT FILAMENT 


The potential of an infinitely long straight filament of constant 
linear density P’ may be found by means of the well-known theorem 
of Gauss, viz., that the integrated value over any closed surface of 
the component of attraction normal to this surface is 4m times the 
mass within. Take as the surface a cylinder of unit length with axis 
coinciding with the attracting filament. The attraction, F, is evi- 
dently everywhere in a direction perpendicular to the filament and 
is a function only of r, the distance away from the filament. Hence 
the surface integral of the attraction over the cylinder of unit length 
is 2rrF and, by the theorem of Gauss, 


2rrF = 4rP’ 
or F=2P"/r (1) 
Then the potential is given by 
¢ = —SFdr = —P'logr? (2) 


where the arbitrary constant of integration is of no consequence in 
our problem since only derivatives of ¢ are involved. 


APPENDIX C 


ADDITIONAL NUMERICAL EXAMPLES 


To show further the method of using the mathematical analysis 
when solving certain practical problems, and also to indicate the mag- 
nitude of the shearing stresses which may be expected, the following 


examples will be solved. In all these cases assume Poisson’s ratio as 


0.25 and FE = 30 000 000 lb. per sq. in. 

Case I, Crossed Cylinders.—Assume R; = 16.5 in., Re = 14 ine 
P = 25 000 lb. This corresponds to the condition of a freight-car 
wheel on a rail having a radius of head of 14 in. 
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B/A = R,/Rz = 16.5/14 = 1.18 


1 1 1 1 
p+a=%(e +g) = 4 (eq + G) - 9.0080 


From Fig. 2, for B/A = 1.18 we have 


IK = @O.08 
IME 
pee Rak: 
ah 
@1 
= = 0.44 
a 


A 
GZ, a hoes = 0.388 


_ 2(1 — o?) i 
~ E(A+B) ~~ 16 X 108 (A + B) 
21 O47 <2 25.000 : 
C= 1.13 X 10? = (() VG tat 


a 0.388 X 0.276 X 107 = 113 000 lb. 
GZ, —i1Y,) = 0.388 — = 2 
A 9.47 per sq. in. 


=. 92472 10-7 


A 


Maximum shearing stress = 4(,Z, — 1Y,) = 56 500 lb. per sq. in. 
b= ka = 0.90 X 0.276 = 0.248 in. 
Area of contact = rab = 3.1416 X 0.276 X 0.248 = 0.215 sq. in. 
Depth to point of maximum shear =0.44a=0.44 X 0.276=0.121in. 
This stress difference may be checked by using the empirical 
equation 


Z - 235500) PP 23.500 25 000 
42 ~ 18 = (Ri / Re)?! Re ~ (1.18)9 V4)? 
23002 any 

= 016 x BEsieat . per sq. in. 


Case II, Crossed Cylinders.—Assume R, = 40 in., R, = 14 in, 
P = 60000 lb. This may be considered as corresponding to a loco- 
motive driving wheel having a load of 30 000 lb., plus 100 per cent 
allowance for impact (as frequently assumed for short spans), on a 
rail having a radius of head of 14 in. 
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B/A = R,/Rz = 40/14 = 2.86 
Vy a a) = 0.0482 
Bet Ay) AV 40 at) ae 
From Fig. 2, 
k -= 0.500 
AP 
— = 0.430 
a 
@1 
== Onl, 
a 
A 
, ae — 1¥y) = 0.274 
i 1 
Naas on = 1279050 10e" 


16(A + B)108 16 X 0.0482 X 10° 


OMsOSCI0 be ee 


42% X 60 000 
a= 

b = 0.5a = 0.283 in. 
Area of contact = rab = 0.502 sq. in. 


a 0.274 X 0.565 X 107 ~—_ 149 600 Ib. 
Agia 12505 per sq. in. 


iZ, ie ay 


li 
= 
bo 
~J 
TS 

| 


Maximum shearing stress = 59 800 lb. per sq. in. 
Depth to point of maximum shear 2; = 0.317a = 0.179 in. | 
Note that this wheel and loading produce a shearing stress only six 
per cent higher than the freight-car wheel, in spite of the fact that | 
the wheel load is 2.4 times as great. 

Case III, Cylinder on Plane.—Assume, as in Case I, radius of © 
wheel of 16.5 in., load 25 000 lb., but assume the head of the rail is 
flat with the wheel making contact on 2 in. of width of the head. 


Width of contact 2b = 0.000564 VPR 


25 000 
0.000564 wor, xX 16.5 


= 0.256 in. 


I 
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Maximum shearing stress 


. 1355°01P 9 1855. 1195500 
= YZ. = i Dy — = \ 


TP MERMS 5 16.5 
= 18 600 lb. per sq. in. 


Depth to point of maximum shearing stress 


0.256 
21 = 0.786b = 0.786 X Ie Sats 0.100 in. 


Case IV, Cylinder on Plane.—Assume, as in Case II, radius of 
wheel of 40 in., load 60 000 lb., but a flat-head rail, with contact 
over 2 in. of width of rail. Load per inch P = 30 000 lb. 


Width of contact 2b = 0.000564 1/30 000 x 40 


= ().618 in. 
1355 .|P 
Shearing stress = vor Wa = 18 600 lb. per sq. in. 


Depth to maximum shearing stress 


0.618 : 
= 0.786 X Se ai 0.243 in. 

Case V, Cylinder on Plane.—In University of Illinois Bulletin 
No. 162,* page 54, are given the results of tests of a segment of a 
roller for the Rall Type Bascule Bridge. This roller had a radius of 
30 in., the Brinell hardness of the material being 160. Using the 
average load per inch (39 800 lb.) given in Table 9 as causing yielding 
in carbon steel, we have a calculated shearing stress of 


= : y. 1355 \" 800 
aximum shear = 4 X 13855 eS 30 


= 24 700 lb. per sq. in. 

No results of physical tests of the material are given in the bul- 
letin, but it. seems reasonable to expect the shearing yield point of the 
material having a Brinell hardness of 160 to be in the neighborhood 
of this calculated stress. 

Case VI, Cylinder on Plane.—In Bulletin 162, just referred to, 
on page 55 are given the loading conditions for the sole plates of the 
Galveston Causeway Bascule Bridge. These plates had'a radius of 


*“Tests of the Bearing Value of Large Rollers,” Univ. of Ill. Eng. Exp. Sta. Bul. 162, 1927. 
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309 in., were 214 in. thick and subjected to a design load of 112 400 | 
lb. per inch of length. For this case: | 
Width of contact 2b = 0.000586 112 400 X 309 = 3.45 in. 


Maximum shearing stress 


1355 | tiot00 ee 
eee BHO . per sq. in. 


Depth to point of maximum shearing stress 
2, = 0.7866 = 0.786 X 1.725 = 1.355 in. 


Based on the shearing stress calculated for this condition one | 
might judge that the plates would be satisfactory, and yet these | 
plates had to be removed from the bridge because they did not stand | ! 
up. Professor Wilson reports that two pieces of these plates, finished |} 
to a thickness of 2 in., having a length parallel to the axis of the cyl- |) 
inder of 4.1 in., and a radius of 309 in. showed permanent set under | 
an average load of 53 000 lb. per in. For this load the calculated || 
shearing stress is 8 900 lb. per sq. in. The failure of the mathematical | 
methods in this case is not difficult to understand. 

It will be remembered that in the derivation of the mathematical | 
equations for the condition of loading of a cylinder on a plane it was | 
assumed that the ratio of the length of contact (parallel to the axis | 
of the cylinder) to the width was infinite, and also that the thickness 
of the bodies was infinite, or at least very large. It is evident, then, | 
that in using the mathematical methods for calculating the stresses | 
existing in this case they are being applied to conditions which are | 
entirely different from the assumptions on which the mathematical | 
analysis was based. i 

Professor Wilson found that, for a block 10 in. deep, length of | 
contact 4 in., radius 309 in., the load required (interpolated from | 
curve) to produce set in a material having a Brinell hardness of 117 
would be about 113 000 lb. per in. In this case also the length of 
contact is small in comparison with its width and hence the calculated 
shearing stress of 13 000 lb. per sq. in. cannot be expected to be 
correct. 
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